Observability of the arrival time distribution using spin-rotator as a quantum clock 
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An experimentally realizable scheme is formulated which can test any postulated quantum me- 
chanical approach for calculating the arrival time distribution. This is specifically illustrated by 
using the modulus of the probability current density for calculating the arrival time distribution of 
spin-1/2 neutral partides at the exit point of a spin-rotator(SR) which contains a constant mag- 
nètic field. Such a calculated time distribution is then used for evaluating the distribution of spin 
orientations along different directions for these partides emerging from the SR. Based on this, the 
result of spin measurement along any arbitrary direction for such an ensemble is predicted. 
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Introduction. — Of late, the question of calculating the 
arrival or trànsit time distribution in quantum mechanics 
has been a tòpic of much interest. For comprehensive 
reviews see, for example, Muga and Leavens and 
Muga et al. 0. A number of schemes Q have been 
suggested in the literaturc for calculating the arrival 
time distribution such as those based on axiomàtic 
approaches, trajectory models of quantum mechanics, 
attempts to define and calculate the arrival time dis- 
tribution using the consistent histories approach, and 
attempts of constructing the time of arrival operator, 
ctc. Thus there is an inherent nonuniqueness within the 
formalism of quantum mechanics for calculating time 
distributions such as the arrival time distribution. 

Against the backdrop of such studies, it remains 
an open question as to what extent these different 
quantum mechanical approaches for calculating the 
time distributions can be empirically discriminated. An 
effort along this direction was made by Damborenea 
et al. who considered the measurement of arrival 
time by the emission of a first photon from a two-level 
system moving into a lascr-illuminatcd region. They 
had evaluated the probability for this emission of the 
first photon by using the quantum jump approach. The 
suitable approximations under which such calculated 
results could be related to Kijowski's axiomàtic arrival 
time distribution and the arrival time distribution 
defined in terms of the probability current density (not 
its modulus) were also discussed. Subsequently, further 
work was done along this direction by Hegcrfcldt et 
al. who made more precise the connection of this 
approach with Kijowski's distribution. 

In this paper we address this question from a new 
perspective so that one can start from an axiomatically 
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defined time distribution and then directly relate it to 
the actually testable results. In order to illustrate this 
approach, here in particular, we use a time distribution 
postulated in terms of the (normalised) modulus of the 
probability current density. Based on this, we derive a 
distribution of spin orientations along different direc- 
tions for the spin-1/2 neutral partides emerging from a 
spin-rotator (SR) which contains a constant magnètic 
field. Such a calculated distribution function can then 
be tested by suitably using a Stern-Gerlach(SG) device, 
as explained later. Thus the scheme formulated in 
this paper can also be viewed as a verification of the 
observability of the quantum probability current density. 

Unlike the position probability density, the status of 
probability current density in quantum mechanics as an 
observable quantity has remained a problemàtic issue; 
see, for example, Kan and Griffin || who pointed out 
that for a many-particle system, any linear operator 
representation for velocity is inconsistent with a linear 
operator representation for the probability current 
density, such as the one constructed by Landau [Ü for 
the quantum thcory of superfluid helium. Nevertheless, 
in the context of singlc-particlc dynamics, the proba- 
bility current density has been used in the quantum 
mechanical predictions of time distributions such as the 
arrival time 0, 0, llfjj . tunneling and reflection times 

m 

Next, let us consider the analysis of the experimental 
tcchniques for measuring the "arrival time" or "time 
of flight" . We note that such analysis is usually done 
semi-classically or classically 01 ■ Therefore it is curious 
that the question of a consistent quantum mechanical 
treatment of the measurement of time has remained 
murky ever since Pauli's argument fl3l | that "there 
cannot be a self adjoint time operator conjugate to 
any Hamiltonian bounded from below" . Subsequently, 
a number of authors |14j have pointed out various 
conceptual and mathematieal problemàtic aspeets of 
this question. On the other hand, several specific toy 
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models |15| have also been proposed to investigate the 
feasibility of how actually the measurement of a time 
distribution can be performed in a way consistent with 
the bàsic principies of quantum mechanics. 

Now, since such debatés arise essentially if one consid- 
ers how to directly measure time in quantum mechanics, 
here we bypass this vexed issue by adopting the following 
strategy. We consider the SR as a "quantum clock" 
where the bàsic quantity which determines the actually 
observable results is the probability density function 
IT (0) which corresponds to the probability distribution 
of spin orientations along different directions for the 
partides emerging from the SR, <p being the angle 
by which the spin orientation of a spin-1/2 neutral 
particle (say, a neutron) is rotated from its initial 
spin polarised direction. Notc that this angle is 
determined by the trànsit time (£) within the SR. Hencc 
the probability density function II (</>) stems from II (t) 
which represents the distribution of times over which 
the partides interact with the constant magnètic field 
whilc passing through the SR. It is the evaluation of this 
quantity II (í) which critically depends on what quantum 
mechanical approach one adopts for calculating such a 
time distribution. 

The plan of this paper is as follows. We will first 
elaborate on the relevant setup, with a discussion of 
how the estimation of the quantity II (4>) can actually 
be tested by using a SG device. Subsequently, we 
will outline the key ingredients of the specific schcmc 
adopted in this paper for calculating II (0) in terms 
of FI (t) using the modulus of the probability current 
density. This will be followed by illustrative numerical 
estimates. 

The setup. — Traditionally, a SR has been mainly used 
for the neutron interferometric studies [TÜI ]. Application 
of the Larmor precession of spin in a magnètic field has 
earlier been discussed, for example, in the context of the 
scattering of a plane wave from a potential barrier |l7j . 
On the other hand, the scheme proposed in this paper 
explores an application of Larmor precession such that 
one can empirically test any given quantum mechanical 
formulation for calculating the arrival time distribution. 

We consider an ensemble of spin 1/2 neutral partides, 
say, neutrons having magnètic moment fx. The spatial 
part of the total wave function is represented by a 
localised narrow Gaussian wave packet ip (x,t = 0) (for 
simplicity, it is considered to be one dimensional) which 
is peaked at x — at t = and moves with the group 
velocity u. Thus the initial total wave function is given 
by * = V (x, t = 0) ® x (* = 0) where x (* = °) is the 
initial spin state which is taken to be same for all 
members of the ensemble. 

The SR used in our setup (Fig. 1) has within it a 



constant magnètic field B = Bz directed along the +z 
- axis, confined between x = and x = d. Within 
the SR, the spatial part of the total wave function is 
assumed to propagate freely, while its spin part interacts 
with the constant magnètic field. This assumption 
is justified in our setup because, for our choices of 
parameters, the magnitude of the Zeeman potential 
cnergy of the interaction of spin of the neutron with 
the constant magnètic field (~ 0.01 neV) is exceedingly 
small compared to the kinetic energy of the neutrons 
(~ 0.01 eV). Hence to a very high degree of aceuracy we 
can consider the evolution of the spatial wave function 
within the SR to be free. Therefore the spatial and the 
spin parts of the total wave function can be considered 
to evolve independent of each other in a tensor product 
Hilbcrt space H = Hí <g> H2 where HI and H 2 are the 
disjoint Hilbcrt spaces corresponding to the spatial and 
the spin parts of the total wave function respectively |l6j | . 

Here in our analysis we assume that the spin part of 
the total wave function begins to interact with the con- 
stant magnètic field in the SR at the instant (t = 0) when 
the peak of the incoming wave packet is at the entry 
point (x = 0) of the SR. Thus the calculational proce- 
dure adopted here is esssentially vàlid for a sufficiently 
narrow wave packet. Then it can be assumed that the 
entire ensemble of partides corresponding to the initial 
wave packet start interacting with the magnètic field at 
t = 0. This assumption is crucial in this schcmc in order 
to enable the arrival / trànsit time distribution l·l (t) to 
be mapped onto l·l (</>). 

It is important to mention here that in the Standard 
approach, the question as regards the initial instant at 
which a propagating wave packet of any arbitrary width 
starts interacting with a localizcd potential (e.g. the lo- 
calized magnètic field in our setup) is intrinsically prob- 
lemàtic since there is no unique criterion for fixing this 
instant. Implications of this nonuniqueness have hithcrto 
remained unaddrcssed in the literature. Following the 
usual procedure (expected to be vàlid for narrow wave 
packets) we have assumed in our paper that the spin 
part of the wave function begins to evolve under a given 
localized potential essentially from the instant when the 
peak of the wave packet reaches the entry point of that 
potential. This criterion is of course not the only crite- 
rion one can use to define this initial instant and one can 
ccrtainly use any other criterion to fix this initial instant. 

Basically there are two different kinds of non- 
uniqueness in the quantum mechanical treatment of our 
problcm; one is that there is no unique criterion availablc 
within the Standard framcwork of quantum mechanics for 
fixing the initial instant when the spin part of the wave 
function starts to interact with the localized potential; 
the other non-uniqueness is inherent within the formal- 
ism of quantum mechanics is regarding the time- dur ation 
over which the wave packet interacts with the localizcd 
potential which in our paper is fixed by the arrival/transit 
time distribution. These are the crucial conditions rclc- 
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FIG. 1: Spin-1/2 partides, say, neutrons with initial spin orienta- 
tions polarised along the +x - axis and associated with a localized 
Gaussian wave packet (peaked at x = 0, t = 0) pass through a spin- 
rotator (SR) containing a constant magnètic field B dirccted along 
the +z - axis. The partides emerging from the SR have a distribu- 
tion of their spins oriented along diflerent directions. Calculation 
of this distribution function is expcrimentally tested by measuring 
the spin observable along a direction n (9) in the xy-plane making 
an angle 9 with the initial spin polarised along +x- axis. This is 
donc by suitably oricnting the direction n (9) of the inhomogeneous 
magnètic field in the Stern-Gerlach (SG — n) device. 



vant to our work. 

The important point is that using our scheme one can 
test experimentally any postulated quantum mechani- 
cal approach for calculating the arrival time distribution 
with diffcrcnt criteria for fixing the initial time of interac- 
tion. Notc that a number of approaches have been sug- 
gested in the literature for calculating the arrival time 
distribution. We have adopted in our paper one par- 
ticular approach, viz., the current density approach to 
calculate the arrival time distribution using a particu- 
lar method of calculation ( e.g. fixing the initial time 
of interaction in terms of the peak), as an example, to 
illustrate our general scheme. 

Next, we recali that when a spin-polarised particle 
(say, a neutron) passes through the constant magnètic 
field within a SR, its spin orientation is rotated by 
an angle <j> with respect to the initial spin polarised 
direction along +x axis. This angle is fixed by the time 
(i) spent by the particle within the SR, given by the 
well known quantum mechanieal relation </> = 2u>t where 

uj = fj,B/n Eg. 

Now, lct us consider an ensemble of partides, say, neu- 
trons passing through the SR where initially all members 
of this ensemble have their spins polarised along, say, the 
+x axis. Given the same initial spin state they evolve 
over different times (characterised by II (£) , the distribu- 
tion of trànsit times within the SR) under the interaction 
with the constant magnètic field within the SR for many 
repetitions of the experiment. 

Here we would like to stress that in our setup, the very 
fact that a distribution of spins emerges from the SR 
implies the existence of a distribution of trànsit times 
within the SR, since the spin rotation is proportional to 
the trànsit time. Hcncc in order that a distribution of 
spins emerges from the SR, the spin part of the wave 



function going through the SR must necessarily interact 
with the SR magnètic field over different times in each 
run of the experiment. It then follows that the unitary 
evolution operators U = exp(—iHt/H) are different 
for each repetition of the experiment, although the 
spin interaction Hamiltonian H is the same for all of 
them. Thus the emergent spin states get polarised along 
different directions and consequcntly the final ensemble 
of partides emerging from the SR is in a mixed state 
of spin states polarised along various directions with 
different respective probabilities. 

Hence we can write the final density matrix of the total 
ensemble at any time which is large enough so that by 
which all the partides of the ensemble (total wave packet) 
have passed through the spin rotator (SR) to be given by 

Wf=Y,m)\x{t)){x{t)\ (i) 
t 

where |x(í))'s oceurring in the right hand side of Eq.(l) 
are the time evolved pure states which have evolved un- 
der the given potential within SR over different times 
(denoted by the symbol "t"). The final density matrix 
by combining these time evolved spin states is written at 
a sufficiently large time by which all the members of the 
ensemble have passed through the SR. Now, using the 
relation t = 4>/2lú, one can rewrite the density matrix 
given by Eq.(l) in the following form 

w f = J2n(<t>)\x(<t>))(x(<l>)\ (2) 

where the summation in Eq. (2) is over the different 
vàlues of (f> corresponding to different vàlues of trànsit 
time (i) within the SR, and \x (<fi)) is the normalised spin 
state which represents the spin polarization along any 
direction making an angle </> with the +x axis. Here II (<f) 
is the normalised probability density of spin orientations 
which is obtained from II (í) through the relation t = 
4>/2oj. The quantity H((f>)d(j) represents the probability 
of spins emerging from the SR having their orientations 
within the angles 4> and <fi + d<fi. Note that using Eq.(2), 
it follows that 

Tr{W f } = ^( Ml |^ / k)=^n(0)^|(u ï | x (0))| 2 

i (f> i 

= £n(0) = i (3) 

where we've introduced the states as a complete 

set of orthonormal basis for any spin state |x (</>)), and 
hence J^i \ ( u í I X (</>)) I = !• Note that the above result 
is vàlid regardless of whether the spin states \x (4')} f° r 
the different vàlues of 4> are orthogonal or not. 

Tcstability of II (</>) using the Stern-Gerlach device. — 
Now, for testing the scheme we have outlincd for calculat- 
ing the probability density function II (</>) , let us consider 
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the measurement of a spin variable, say ag, by a SG dc- 
vice (Fig.l) in which the inhomogeneous magnètic field 
is oriented along a direction n (9) in the xy-plane making 
an angle 9 with the initial spin-polarised direction (+x 
axis) of the partides. Then for the spins of the partides 
emerging from the SR polarised along different directions 
(with the probabilities II (0)) making angles with the 
+x axis, the probabilities of finding the spin component 
along +9 direction and that along its opposite direction 
are respectively given by 



p+V) 



P_(6) 



n (0) Cos 2 



n (0) Sin 



2(0 



(4) 



(5) 



where P+ (9) + P- (9) = 1, and here we are essentially 
restricting to the situations in which the relevant 
parameters d, u and B are such that the spin rotation 
angles f° r a U the partides emerging from the SR 
are restricted between = and = 2ir. To cxplain 
in morc dctail how Eqs.(4) and (5) are derived, let 
us first consider partides passing through the spin 
rotator (SR) with all their spins oriented along a 
definite direction, say, x-axis. The initial x-polarised 
spin state can be written in terms of the z-bases | f ) 
and | | ), as X (0) = 1/V2 (| T ) z + \l >,)■ Then in 
such a case, the spin polarised state rotates only in 
the xy-plane. If is the rotation angle with respect 
to the initial spin orientation along x-axis, such a 
rotated spin state in the xy-plane can be typically 
written as x{<t>) = VV2 (| T }, + e^| I )J. Now, for 
the purpose of measurement after the spins emerge 
from SR, if onc applies SG-magnetic field along a 
direction n (9) in the xy-plane making an angle 9 with 
the x-axis, then the bases states for the spin operator 
Ò e are respectively | |)e = (| T ) z + e t6 \ j )J 

and | l) = 1/V2 (| T) 2 -e l9 | l) z ). Then for this 
spin measurement the probabilities of getting | D 6 and 
| l)o arc P+ (9) = | e (t | x(0)>| 2 = Cos 2 {9 - <f>)/2 and 
P-(8) = \e(i I X(4>))\ 2 = Sin 2 (9-cb)/2 respectively 
Now, since in our setup, instead of a definite spin 
polarised state, we have considered a distribution of spin 
orientations along different directions for the partides 
emerging from the SR characterised by the distribution 
function 11(0), conscqucntly we have the expressions for 
the probabilities P+{9) and P-{9) as given by Eq.(4) 
and Eq.(5) respectively. Notc that the departures 
obtaincd from the prediction given by semiclassical 
approach (which assumes that, neglecting the effect due 
to the spreading of the wave packet, the spins of all 
members of the ensemble rotate by an amount which is 
determined by the time spend by the peak of the wave 
packet in traversing the region within the SR) can be 
observed by sensitive measurements of P+ (9) and P_ (9). 



It is these probabilities P + (9) and P_ (9) which 
constitute the bàsic observable quantities in this scheme 
which are determined by the distribution of spins II (0) 
of the partides emerging from the SR. The estimations 
of these probabilities crucially depend on how one 
calculates the quantity II (0) whose evaluation, in turn, 
is contingent on the procedure adopted for calculating 
the relevant time distribution II (i). As mentioned 
carlicr, the specification of such a time distribution 
is not unique in quantum mechanics. For the setup 
indicated in Fig. 1 , II (i) represents the arrival time 
distribution at the exit point (x = d) of the SR, which 
is also the distribution of trànsit Urnes (t) within 
the SR. In the specific scheme we are using, II (í) is 
taken to be representcd by the modulus of the prob- 
ability current density |J (X, í)| (suitably normaliscd) 
evaluated at the spatial point X(x=d,y=0,z=0); i.e., 
we take II(í) = |J(X, t)\/ J °° |J(X, t)\dt. Leavens 
0, has justified the above interpretation of the 
modulus of probability current density |J(X, í)| as 
an (unnormalizcd) arrival time distribution using the 
Bohm's causal model of quantum mechanics. Hencc 
n(0) = |J(X, 0)|/ / 2?r |J(X,0)|#. Thus in this scheme, 
the calculation of 11(0) ultimatcly hinges on cvaluating 
J(X,0) from J(X,í). 

The evaluation of 11(0). — Since we have to first cal- 
culate J(X,í), we begin by rccalling that the Standard 
expression for the non-relativistic quantum probability 
current density is given by 



J(x,í) = Re 



r (x,t) 



ih 



VV> (x,í) 



(G) 



which satisfies the quantum mechanical cquation of con- 
tinuity given by 

| + V.J = (7) 
where the position probability density p (x, í) = 

r (x,í)v(x, t). 

Then comes a key point. If one adds any divergence- 
free term to the above expression for J(x, í), then 
the new expression also satisfies the same equation of 
continuity. Hence there is a nonuniqueness inherent in 
the nonrelativistic expression for the probability current 
density. Curiously, this point has not been noted even in 
the premier textbooks like that by Landau and Lifshitz 
[ï^ . Merzbacher J2£j. However, relatively recently this 
problem of nonuniqueness has been highlighted [2l| and 
it has been pointcd out that the probability current 
density derived from the Dirac equation for any spin- 1/2 
particle is unique and even in the non-relativistic limit 
it contains a spin- dependent term which is present in 
addition to the expression for J(x,í) given by Eq.(6). 
Interestingly, one can further arguc that this property 
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of the uniqueness of probability current density is not 
specific to the Dirac equation, but is a consequcncc 
of any relativistic quantum mcchanical equation. The 
argument is as follows. 

The probability current density obtained from any 
consistent relativistic quantum mcchanical equation 
needs to satisfy a covariant form of the continuity 
equation of j M where the zeroth component of is 
associated with the position probability density. If one 
replaces j M by j which is also conserved, i.e., d^j = 
where = + a M (a M is an arbitrary 4- vector), then 

the zeroth component of j 1 * Q ) will have to be the same 
as the position probability density given by j°. Hence it 
follows that a° = 0. 



Next, we consider this current as seen from another 



Lorentz frame. This is given by j 



-0/ 



Hence 



in this frame j = j + a , and again if the position 
probability density has to remain unchanged, then one 
must have a ' = 0. But wc know that the only 4-vector 
whose fourth component vanishes in all frames is the 
null vector. Thus a M = 0. It thereforc follows that for 
any consistent relativistic quantum mcchanical equation 
satisfying the covariant form of the continuity equation, 
the relativistic current is uniquely fixed. This uniqueness 
is also preserved in the non-relativistic limit of the 
relevant relativistic equation. 

Now, in order to makc this paper self-containcd, we 
briefly recapitulate the kcy steps involved in deriving the 
expression for the probability current density in the non- 
relativistic limit from the Dirac equation in 3 + 1 dimen- 
sion for a free spin-1/2 particle of rest mass tuq given 
by 



9* 



hc 9 2 
— oíí- \~ pm c 

l OXs 



* 



(8) 



where 



^ 
o-i 



-I 



*i 

* 2 



and í'l, *2 arc individually two component spinors. Sub- 
scqucntly, Eq. (8) lcads to two coupled cquations 



9*i 



9*j 
dt 



9* 2 

dxi 



dxi 



-m c^i 
n 



T m c * 2 
n 



(9) 



(10) 



Then taking the positive energy solution *2 oc 
exp(—iEt/H), substituting it in Eq.(10) and putting 
E = m c 2 in the non-relativistic regime, we get 



*; 



ih .9*1 
-cr, 



2moc dxi 



(11) 



Multiplying Eq. (9) by *| from the left and multiplying 
again the hermitian conjugate of Eq. (9) by *i from the 
right, we add the two equations. Substituting the value 
of *2 from Eq.(ll) in this resulting equation, one can 
then obtain the following equation given by 



|(*rt*i 



(12) 



d 

dxi 
= 



ih 

2mn 



*iVí 



9*i 

1 dxi 



'9*i+ 
dxi 



cr,*! 



Now, comparing Eq.(12) with Eq.(7), it is seen that the 
Dirac current in 3+1 dimension for a free spin-1/2 neutral 
particle in the non-relativistic limit is of the form given 
by 



J (x, t) 

ih 



(13) 



2toq 



*Íctí 



9*i 



9*J 



cr,*! 



where *i is a two component spinor which can be written 
as *i = ip (x, í) x (í) ■ Simplifying Eq.(13) by using the 
Gordon decomposition |2Ïl l22j . one finally obtains 



J(x,í) 

= Re 

= 3sch ( x , t) + Jspin ( x , t) 



m 



m 



(14) 
[Vp x s (í)] 



where p = \ip (x, í)| 2 , s(í) = f X^(*) <t x( í )j a = cr x x + 
crj,y + cr z z and x(í) x{t) ~ 1- I n Eq.(14), the first term 
represents the usual Schroedinger current (J5 c /,(x,í)) 
and the second term gives the contribution of the spin 
dependent current (Jspinfat))- The above decomposi- 
tion is possible because there is no spatial dependence on 
the spin state x(t) which is only time dependent here in 
our case. 

Next, in the context of our setup, in order to evaluatc 
the quantity J (x = X, <fi) , we first consider the spatial 
part of the total wave function. As mentioned earlier, it 
is taken to be a one dimensionat Gaussian wave packet 
which is peaked at the entry point (x = 0) of the SR at 
t = (Fig. 1). The initial spatial wave function is then 
given by 



ip(x,t = 0) = 



■ exp 



ikx 



(15) 



where cto is the initial width of the associated wave 
packet. The wave number k = rriQu/K where u is the 
group velocity of the wave packet moving along the +x 
- axis. Since the spatial part of the total wave function 
propagates freely, being unaffectcd by the constant mag- 
nètic ficld confined within the SR, the Schroedinger time 
evolved spatial wave function calculated from the initial 
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wave function tjj(x,t = 0) given by Eq. (15) is of the 
form 



ifi(x,t) 



1 



(2ttA?) 



1/4 



exp 



utf 



4A t a 



ik(x 



1 



-ut) 



(16) 

and hence the time evolvcd position probability dcnsity 
is given by 



p(x,t) 



1 



exp 



(x — ut) 



(17) 



where A t = o~o ( 1 

,2.2 \l/2 



//// 



2mocg 



and o t = \AÀ = 



(T 1 



; o't is the width of the wave packet at 



any instant t. 



Next, we consider the spin part of the total wave func- 
tion. As mentioned earlier, the initial spin of a spin-1/2 
neutral particle is taken to be polarized along the +x- 
axis; i.e., the initial spin state is given by 



X(0) 



1 



V2 



[|í>, + U>, 



and hence 



s(0H-X f (0)<r X (0) 



(18) 



(19) 



Then we proceed to calculate the time evolved spin part 
of the total wave function undcr the interaction Hamil- 
tonian H = /íct.B. For this purpose we note that the 
constant magnètic field B = Bz within the SR, confincd 
between x = and x = d, is directed along the +z— axis. 
Then the time evolved spin state x (t) is given by 



x(t) = 

whcncc 



exp 



-il·lt 



X(0) 



1 



"li),] 

(20) 



s (í) = 2 X f (í) VX (t) = ^ ( cos 2ujt x + sin 2wí y) (21) 

where w = pB/h. Given the above expressions for 
tp(x,t), p(x,t) and s(t) corresponding to Eqs. (16), 
(17) and (21) respectively, one can now calculate the to- 
tal probability current density at any given point from 
Eq.(14). We specifically evaluate it at the exit point 
X(x = d, y = 0, z = 0) of the SR. Then the expres- 
sion for the total current density J(X,í) reduces to the 
form given by 



J(X,t) 

, . f (d - ut) h 2 t } _ 

. , . f h (ut — d) . 

+ p(x = a, í) < ^— sm 2cjí > z 

[ 2mo(Tí 



(22) 



The first term in Eq.(22) is the usual Schroedinger cur- 
rent, and the second term represents the additional con- 
tribution arising from the spin of the particle. Next, 
substituting t = 4>/2lü in Eq.(22), one gets the following 
expression for the probability distribution of spin orien- 
tations for the partides emerging from the SR given by 



J(x,0) = P (d,4>)lu + 



uó 



2u 



4 2 4 , Er. 



p (d, </>) 



■{%-*) 
2m al 



sm i 



Jsch (d, 4>) 



where p (d, (fi) 



■ exp 



J Spin {d, (fi) 



(23) 
and U(f) = 



TO 



1 



. Thus the probability density func- 



tion corresponding to the distribution of spins emerging 
from the SR is given by 



n(0) 



\3(X,<fi)\ 



C |J (X, 0)|. 



(24) 



where |J (X, cfi) \ is calculated from Eq.(23). Next, we pro- 
ceed to present the results of a few numerical estimates 
for the observable probabilitics P+ (9) and P_ (8) deter- 
mined by Eqs. (4) and (5) respectively, where II ((fi) is 
calculated using Eqs. (23) and (24). 

Numerical estimates for P+(0) and P_(0). — These 
estimates have been done using different vàlues of the 
relevant parameters; viz. the initial width (oo), the 
spatial extension (d) of the constant magnètic field 
confincd within the SR, the group velocity (u) of the 
peak of the wave packet and the magnitude of the 
constant magnètic field (B). 

The choices of these parameters in our calculations 
are constrained by the condition that the parameter d 
should be sufficiently large compared to the half width 
of the wave packet peaked at the exit point (x = d) of 
the SR; also, the relevant parameters d , u and B are 
chosen such that the spin rotation angle (fi is effectively 
confincd between = and (fi = ir so that the SG 
magnet does not come in the way of the neutron beam. 

Now, for presenting the results of our estimates, we 
first show a few representative curves (Figs. 2a and 2b) 
which correspond to the probability density functions 
II (efi) of the spin orientations of the partides emerging 
from the SR. Note that the curves in Figs. 2a and 2b are 
respectively associated with two different sets of choices 
(say, I and II) of the parameters d, u, and B. Whilc 
the two curves in Fig. 2a correspond to two different 
vàlues of (T (10~ 5 and 10 -4 ), the two curves in Fig. 2b 
correspond to those vàlues of ao . The curves in Fig. 2a 
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FIG. 2: The quantity Tl(<f>) denotes the probability density func- 
tion which represents the distribution of spin orientations of the 
partides emerging from the SR. This quantity is plottcd against the 
angle of spin rotation cf>. The two curves in Fig.2a correspond to the 
set I of the choices d = lcm, n = 3x 10 5 cm/s and B = 10 gauss; 
both these curves are peaked at <f> = <f>i = 34.94767°. The two 
curves in Fig.2b correspond to the set II of the choices d = 2cm, 
u = 3 X 10 5 cm/s and B = 10 gauss; both these curves are peaked 
at (f> = 02 = 69.89534°. 



for the set I correspond to d = lcm, u = 3 x 10 5 cm/s 
and B = 10 gauss, while the curves in Fig. 2b for the 
set II correspond to d = 2cm, u = 3 x 10 5 cm/s and 
B = 10 gauss. 

Note that both the curves in Fig. 2a represent the 
probability density function II (</>) peaked at cf> = (f>i = 
34.94767°, while the curves in Fig. 2b represent II (0) 
peaked at (j> = cf>2 — 69.8953°. For different choices of 
the initial width (oo) with the other relevant parameters 
d , u and B remaining fixed, it is seen from each of Figs. 
2a and 2b that the qualitative nature of these curves and 
the location of their peak remain the same, while their 
varianecs differ, increasing with decreasing vàlues of <jq. 

Next, a fcw rcprcscntativc results of the numcrical 
computations based on using Eqs. (4), (5) and (24) are 
given in Tables 1 and 2 which correspond respectively to 
the sets of vàlues I and II of the parameters d ,u and B 
(corresponding to the Fig. 2a and Fig. 2b respectively). 
The results shown in any one of these Tables indicate 
how the observable quantities P+ (9) and P- (9) vary 
with different initial widths <7q, corresponding to a 
specific orientation (9) of the inhomogencous magnètic 
fleld in the SG devicc (Fig.l). 

It is seen that for a given value of 9, for both the sets I 
and II of the choices of the relevant parameters d ,u and 
B, the variation in the vàlues of the probabilities P+(6) 
and P- (9) is very small, but is detectable for smaller 
vàlues of <tq. Morc comprehensive estimates of P + (9), 
P- (9) and fluctuations in the distribution of spins emerg- 



ing from the SR, for a wider variation of 9 and other rel- 
evant parameters d ,u and B, will be presented in a later 
study. 

Table. 1. The quantities P + (0)and P_ (0)denote proba- 
bilities of the spin measurement along different directions n(8) 
making angles 8 with the initial spin polariscd along +x - axis. 
The numerical vàlues of (9) and P_ (d) are calculated for 
different initial widths cro of the Gaussian wave packet. In this 
Table, the results are presented for three different vàlues of d for 
the set I of the vàlues of the other relevant parameters d = lcm, 
u = 3 X W s cm/s and B = 10 gauss, while 0i = 34.94767° at 
which the curve of IT (0) is peaked. 





8 = cf>i 
= 34.94767° 


t 


) 

+ 60° 


( 


) 

+ 90° 


^0 

(cm) 


p+V) 


P-{0) 


p+(o) 


P_(6) 


P+(0) 


P-(0) 


10" b 


1.00000 


0.00000 


0.75000 


0.25000 


0.50000 


0.50000 


10~ tí 


1.00000 


0.00000 


0.75000 


0.25000 


0.50000 


0.50000 


io- v 


0.99998 


0.00002 


0.75002 


0.24998 


0.50003 


0.49997 


io-« 


0.99886 


0.00114 


0.75242 


0.24758 


0.50345 


0.49655 



Table. 2. The quantities P + (6>)and P_ (d) denote proba- 
bilities of the spin measurement along different directions n(6) 
making angles 9 with the initial spin polarised along +x - axis. 
The numcrical vàlues of P_|_ (6) and P_ (6) are calculated for 
different initial widths <ro of the Gaussian wave packet. In this 
Table, the results are presented for three different vàlues of for 
the set II of the vàlues of the other relevant parameters d = 2cm, 
u = 3x W 5 cm/s and B = 10 gauss, while <j> 2 = 69.89534° at which 
the curve of IT (0) is peaked. 





8 = <t>2 

= 69.89534° 


9 

= ^2 + 60° 


9 

= (t> 2 + 90° 


(cm) 


P+(0) 


P-{0) 


P+{0) 


P-(0) 


p+V) 


P-(0) 


10" b 


1.00000 


0.00000 


0.75000 


0.25000 


0.50000 


0.50000 


10~ tí 


1.00000 


0.00000 


0.75000 


0.25000 


0.50000 


0.50000 


io- v 


0.99995 


0.00005 


0.75004 


0.24996 


0.50006 


0.49994 


io-« 


0.99546 


0.00454 


0.75355 


0.24645 


0.50672 


0.49328 



Summary and Outlook. — In the setup discussed in 
this paper, the observable quantities are the probabilities 
P+ (9) and P- (9) which correspond to the measurement 
of a spin variable along any direction by a SG device 
performed on the partides emerging from the SR. 
Evaluations of these quantities crucially depend on the 
probability distribution II (0) of the orientations of spins 
of the partides emerging from the SR. This in turn 
depends on the quantity 14 (t) which corresponds to the 
distribution of trànsit times over which the partides 
interact with the magnètic field while passing through 
the SR. 

The quantity II (t) is calculated in this paper for spin- 
1/2 partides in terms of the modulus of the probability 
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current dcnsity. Hcnce the estimates presented here 
for the observable probabilities P + (8) and P_ (8) are 
ultimately determined by the modulus of the probability 
current density. Thus if the experimental results for 
such a setup corroborate such predictions, this would 
constitute a verification of the observability of the 
probability current density. 

As mentioned carlicr. bcacuse of an inherent 
nonuniqueness, there are also other quantum me- 
chanical approaches which can be used to evaluate the 
time distribution II (t) , apart from the specific schcmc 
we've used in this paper based on the modulus of the 
probability current density. In the context of our setup, 
it should be instructive to derive the respective pre- 
dictions for the observable probabilities P + (8), P- (8) 
from thcse diffcrcnt approaches with different criteria for 
fixing the initial time at which the spin part of the wave 
function starts interacting with the localized magnètic 
field. One can then compare those predictions with the 
results of the actual experiment. Thercforc it seems that 
an experimental study of the cxamplc discussed in this 
paper could be a worthwhilc effort. 

An interesting question of both theoretical and of ex- 
perimental relevance is as to what effect on the analysis 
does the assumption of a Gaussian wave packet have. 



It should definitely be instructive to study the situation 
taking a non-Gaussian wave function for the spatial part 
which is proposed to be done as a sequel to this work, 
apart from other sequels based on different criteria for fix- 
ing the initial time of interaction and the duration over 
which a propagating wave packet interacts with the lo- 
calised potcntial. 
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